Detecting and localizing the foci in human epileptic seizures 
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We consider the electrical signals recorded from a subdural array of electrodes placed on the pial 
surface of the brain for chronic evaluation of epileptic patients before surgical resection. A simple 
and computationally fast method to analyze the interictal phase synchrony between such electrodes 
is introduced and developed with the aim of detecting and localizing the foci of the epileptic seizures. 
We evaluate the method by comparing the results of surgery to the localization predicted here. We 
find an indication of good correspondence between the success or failure in the surgery and the 
agreement between our identification and the regions actually operated on. 



Epilepsy is defined as a condition of chronic un- 
provoked seizures. It is one of the most common 
neurologic disorders, affecting more than 50 mil- 
lion people worldwide. About 1/3 of cases are 
not adequately controlled by medication, caus- 
ing a severe reduction in independence and qual- 
ity of life. If an individual's seizures are stereo- 
typed and appear to have a focal onset, surgical 
resection may be a consideration. Localizing ac- 
curately the epileptic focus (or foci) is therefore 
essential part of surgical planning. We propose a 
simple and computationally fast method based on 
phase synchronization, that is able to correctly lo- 
calize the foci of epileptic seizures. The method 
was applied to the data of three patients. The 
overlap of the location of the electrodes identi- 
fied by our method and the resected areas was 
consistent with the degree of surgery success. 



I. INTRODUCTION 

Epilepsy is defined as a condition of chronic unpro- 
voked seizures. It is one of the most common neurologic 
disorders, affecting more than 50 million people world- 
wide. About 1/3 of cases are not adequately controlled 
by medication, causing a severe reduction in indepen- 
dence and quality of life. Seizures often begin in one 
brain region and rapidly spread to other areas of cor- 
tex, manifesting in a wide variety of symptoms, ranging 
from a brief staring spell (absence seizures), to general- 
ized convulsions (grand mal seizures). The mechanisms 
of seizure spread are not yet well understood. If an in- 
dividual's seizures are stereotyped and appear to have a 
focal onset, surgical resection may be a consideration. 

As part of the invasive surgical evaluation, arrays of 
electrodes are chronically implanted onto the surface of 
the brain (FigJTj) in order to record the onset and evolu- 
tion of several spontaneous seizures. The ensemble of sig- 
nals [called the electrocorticogram (ECoG)] is simultane- 
ously recorded and analyzed for detecting the proper cor- 
tical area to be resected [1] . Decisions regarding which ar- 



eas of cortex to resect are made by board-certified epilep- 
tologists, based on careful review of the onset and spread 
of several spontaneous seizures. The goal is to remove 
the epileptogenic zone, the minimal area of cortex that is 
necessary and sufficient to prevent further seizures from 
recurring, if it is not located in eloquent cortical areas 
that are critical for normal function. 

The clinical challenge is how to analyze the ECoG sig- 
nals to achieve a simple, automatic, data-driven, fast and 
reliable identification of the epileptic focus. Tradition- 
ally, linear measures, such as correlation coefficients, were 
utilized to complement the visual inspection. In the last 
two decades methods based on the nonlinear measures 
were introduced to the field. 

It has been shown that a significant reduction of the 
dimension of the system 0, @] (as measured by the corre- 
lation dimension [4| and by the neuronal complexity loss) 
occur in the vicinity of the epileptic focus, not only dur- 
ing the seizures (ictal periods), but also between seizures 
(inter-ictal periods). 

The interdependence between signals from different 
sites, studied by means of nonlinear cross-predictability 
and mutual nonlinear predictability [H, @, also allows 
the location of the focus. Surrogate-corrected nonlinear 
measures [7j were shown to correctly identify the focal 
hemisphere in the inter-ictal recordings. Many of these 
methods involve the reconstruction of the phase space 
by time-delay embedding [8j]. Such measures calculated 
from noisy non-stationary time series usually can not be 
directly related to the dimensionality of the attractor @ 
and the calculations are computationally costly. 

Recently, concepts from synchronization analysis were 
brought to the study of ECoG signals [jj [TO, El E2, ES 
EH El] 5 leading to a series of promising results on local- 
ization and predictability of epileptic seizures. In most of 
these works the increase in the synchronization was found 
only during the seizure, while synchronization clusters, 
identified in [l3j], remain synchronized also during inter- 
ictal periods. Similar stability of local hyper-synchronous 
regions was found also in [14]. In spite of demonstrated 
success, formulation of these methods involve some as- 
sumptions, that we would like to relax in our proposed 



approach, as it will be explained in the next Section. 

The aim of this paper is twofold. First, we describe a 
model-free and fast approach, based on phase synchro- 
nization, that is able to identify the epileptic focus. Sec- 
ond, we apply the method to ECoG data of three pa- 
tients, and show that there is a good correspondence be- 
tween the success or failure in surgery and the overlap 
between the area spanned by the electrodes identified by 
our method and that of the brain tissue actually removed. 



II. METHOD 
A. Patients and ECoG recordings 

In this paper we analyze the ECoG data from three 
female patients with intractable partial complex seizures 
with secondary generalization, referred to as patients A, 
B and C, respectively. 

Patient A is a 17 years old and had her first seizure 
at age six. Her magnetic resonance imaging (MRI) was 
normal, but a positron emission tomography (PET) scan 
revealed reduced metabolism of the left temporal and 
parietal lobes. Spontaneous seizures originated in the 
temporal lobe and spread posteriorly. A total of 80 elec- 
trodes were implanted in 4 arrays, with 16 eletrodes at 
the temporal lobe. A left temporal lobectomy was per- 
formed, leaving her seizure-free for the last two years 
(Engel class I, see Ref. [20| for the classification of the 
surgery outcome). 

Patient B is a 12 years old with a history of develop- 
mental delay and seizures since 8 month of age. Before 
surgery she experienced as many as 30 seizures/day. Al- 
though her MRI was normal, a PET study revealed gen- 
eral hypometabolism of the cerebral cortex. 96 electrodes 
in 5 arrays were implanted, with 32 electrodes in frontal 
grid (F), 32 electrodes in parietal grid (P) and 24 elec- 
trodes at the temporal lobe in three strips (AT, ST and 
PT, cf. Fig. [2]). She received a left frontal and occipital 
topectomy. She had an Engel Class III outcome. 

Patient C is a 7 years old and the outcome was Engel 
class II 16 months after operation, while 3 years after 
surgery she has frequent seizures. She was implanted 
with 96 electrodes, 64 of which were in the parietal grid, 
and 32 electrodes at the temporal lobe. 

For all patients all subdural electrode channels were 
referenced to a scalp electrode on the midline (CPz) and 
then the grand mean was subtracted from each channel 
to remove the contribution of the reference electrode. 

We have analyzed a total of 24 hours of ECoG record- 
ings (mean 8 hours, range 200 minutes- 15 hours 30 min- 
utes), containing a total of 12 seizures (mean 4, range 
3-5). The data were recorded with a sampling frequency 
of 400 Hz and band-pass filtered offline to 0.5-50 Hz us- 
ing moving overlapping windows and a finite impulse re- 
sponse filter with a Kaiser window. 
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FIG. 1: (Color online) Subdural ECoG grid of electrodes 
placed on the surface of the brain for chronic evaluation of 
epileptic patients before surgical resection. 




FIG. 2: (Color online) The computed tomography (CT) im- 
age with marked electrodes (a) and the electrode scheme (b), 
projected on the plane, for patient B. The arrays of electrodes 
are labeled according to their position: F- frontal, P- parietal, 
AT-anterior temporal, ST - sub-temporal, PT - posterior- 
temporal. The same naming convention is used for all pa- 
tients. For description of other marks see discussion after Eq. 

m 
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FIG. 3: Voltage traces of ECoG signals. Three 10 s portions 
of voltage traces of 7 electrodes, taken from a multielectrode 
recording of 96 electrodes of patient B, cf. Fig. [2] The 
seizure onset is marked by a vertical line. The analysis below 
indicates that the signals coming from electrodes F6, F14, 
F15 and P50 belong to the epileptic foci, whereas the F13, 
P29 and AT5 do not. The voltage scale of 1 mV and the time 
scale of 1 second are shown below the figure. 



B. Synchronization Strength Diagrams 

To see the main concepts behind the method, consider 
a one-dimensional oscillatory signal x(t). The associated 
phase space can be fully reconstructed by using the co- 
ordinates x(i), dx(t)/dt, d 2 x(t)/dt 2 , with as many 
derivatives as needed to exhaust the dimensionality of the 
system. In such a reconstructed phase space one can then 
define the Poincare section by noting the points obtained 
when the orbit crosses the surface dx(t)/dt = 0. Obvi- 
ously, this surface coincides with the points of minima or 
maxima in the original one-dimensional signal x(t). In 
this work we focus on the points of maxima, and define 
the phase <j>(t) of the trajectory by increasing its value 
by 2tt after each maximum. Precisely, between the kth 
maximum occurring at t = tk and the (fc + l)th maximum 
occurring at t^+i? we define the phase over the period of 
time tk <t < tfc+i using the linear interpolation [h| G3 

</>(t) = 2nk + 2tt (t k <t< t k+1 ) . (1) 

Note that the resulting monotonically increasing phase 
does not contain any information about the amplitude of 
the original signal and thus the results of the analysis are 
not affected by the amplitude variation [18[. 

Phases can be affected by instrumental artifacts; we 
have carefully ascertained that any signal segments that 
appear identical in all recording channels were manu- 
ally excluded from the signal prior to calculation of the 
phases. It is important to remark that our data feature 
broad band spectra. Alternative measures of the phase 
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FIG. 4: Panels (a),(c) and (e): Synchrograms for 3 neighbor- 
ing pairs of electrodes along [F13-F14 (a), F15-F14(c)] and 
across the grid [F6-F14 (e)], calculated for the reference elec- 
trode F14 (the signals for F6, F13, F14 and F15 are shown 
in Fig(3}. Since the synchronization line appears at ipi — 0, 
a small phase offset of (j) s h = tt/5 was added to the phase 
cj) y for clarity. Panels (b),(d), and (f) show the histograms of 
the phase points distribution within the corresponding syn- 
chrogram. Clearly, only one strong line is present in each 
synchrogram, at the position of the offset phase (tt/5). This 
corresponds to the 1:1 synchronization without visible phase 
difference. 

(such as the one based on the Hilbert transform) im- 
ply the presence of a prominent spectral peak [17]. We 
would like to avoid any a priori assumptions on the spec- 
tral structure of the signals, therefore we opted for the 
above definition of phase, Eq. [lj 

The notion of phase synchronization [h| [TtJ refers to a 
process wherein a locking of the phases of weakly coupled 
oscillators is produced, without implying a substantial 
correlation in the oscillators' amplitudes. The same con- 
cept can be related to our case, for measuring the degree 
of phase synchronization of two different signals, say x{t) 
and y(t). Using the basic definition (pQ), one can indeed 
produce the phases of each signal, denoted as (j) x (t) and 
y (t), respectively. Taking x(t) as the reference signal, 
one can further consider the time points of its maxima tk, 
and read the phase (j) y (tk) at these time points. One then 
defines the following reduced phases ijj m (tk) according to 

^(tk) = <t>y(tk) mod27r , 
fo(tk) = $y(tk) mod47r , 

... , 

ipm(tk) <t>y{tk) mod 27rra . (2) 

The graphic representation of ^m(tk) vs. tk is called a 
synchrogram of order m. We will consider the signals x(t) 
and y(t) as synchronized at ratio n : m when the mth or- 
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FIG. 5: Histograms for synchrograms of orders m = 1, 2 and 
3 calculated for the pair of electrodes F14-F15 (cf. Fig. 2] c 
and d). The dashed lines mark offset positions 7r/5,27r + 7r/5 
and 47T + 7r/5. The same bin size ( 27r/100) was used in all 
histograms. The number of strong lines equals to the order 
of the synchrogram, confirming the 1:1 synchronization. The 
position of the synchronization lines coincide with the value 
of the offset phase within ±7r/50. 
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FIG. 6: The vicinity of the 1:1 synchronization line for syn- 
chrogram of electrodes F14 and F15, with the phase offset 
4>sh = 7r. Note that the synchronization line is clearly iso- 
lated. 



der synchrogram exhibits n distinct horizontal lines (i.e. 
fixed reduced phase ip m (tk) for a series of time Syn- 
chrograms, indeed, have been proved to be useful to vi- 
sualize and trace transitions between different ratios of 
phase locking in biological data, as the cardiorespiratory 
system of a human subject [19[, where non-stationarity 
is often strong. Note that ip m (tk) is, in general, not sym- 
metric, and, within the same pair of signals, may depend 
on which electrode is taken as a reference [17j. An addi- 
tion of any constant phase (/> s h to <j> y (t) does not change 
the synchronization ratio, but merely shifts the position 
of the synchronization lines. Such an offset may facilitate 
the visualization of the line when it appears close to the 
edge of the synchrogram. 



To exemplify this notion in the present context, we use 
the data of the patient B. In Fig [2] we show the CT image 
with the electrode positions marked on it and the plane 
projection of the electrode placement scheme. In Fig. [3] 
the three 10 sec portions of signals from seven electrodes 
are plotted. These portions of the signal correspond to 
the time near the seizure onset. The signals from five of 
these electrodes are used in Figs. HKJ 

We present in Fig. [4] synchrograms constructed from 
the signals associated with three neighboring pairs of 
electrodes, which display the presence of 1:1 phase syn- 
chronization. The fraction of the time, during which the 
signals are synchronized, differs significantly, albeit the 
same geometrical distance between these pairs of elec- 
trodes. 

To test for the presence of more complicated synchro- 
nization ratios, we calculated the synchrograms of higher 
orders. The histograms for them are shown in Fig[5j The 
number of strong lines in each histogram is equal to the 
order of the corresponding synchrogram, thus confirm- 
ing the dominant 1:1 synchronization ratio. Note that 
this ratio is usually assumed in the methods involving 
phase difference [j, [l3|, [H|- We do not make such an 
assumption. In fact, any synchronization ratio may be 
treated on the same footing. In the following analysis we 
concentrate on the first order synchrograms ipi with the 
phase offset s h = it. The phase points that belong to 
the synchronization line may be easily extracted from the 
synchrogram, since they are clearly isolated from other 
points (see Fig[6j). 

The present method shares some similarities with the 
so-called event-synchronization technique [11], though 
this latter strategy is, in general, unable to distinguish 
between different types of n : m lockings. The infor- 
mation contained in the synchrograms can be easily em- 
ployed to extract the level of phase synchronization be- 
tween all the pairs of electrodes at our disposal. To this 
aim we define the statistical strength of synchronization 
Sij(t) (SS) of every given electrode i to any other elec- 
trode j: we take the synchrogram of electrodes i and j, 
first divide the time axis into windows (in our case of 10 
seconds), focusing on the line ip(t k ) = 0. In each such 
window, we calculate the number Nij of phase points 
within the range </> s h =L 0.01 . For i ^ j Sij(t) is defined 
as Nij normalized by the total number of phase points 
throughout the measured time interval, with t being as- 
signed to the first time-point in the window. For i = j 
we set Sij = for all times. Note that this procedure 
does not require time averaging and the choice of the 
time window is dictated by the signal. Precisely, if the 
time window is too large, short events of strong synchro- 
nization will be smeared out and not properly detected, 
thus reducing significantly the sensitivity of the method. 
On the other hand, a too short time window would result 
in a too small number of points over which the fraction 
is calculated, thus introducing strong noise in the eval- 
uation of the synchronization strength. In the present 
case, we have chosen a window length of 10 seconds, that 
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FIG. 7: (Color online) A typical SSD for a focal electrode 
[here, F14 (i=78)]. The electrodes here are numbered and 
correspond to: 1-8 (AT1-AT8), 9-12 (ST1-ST4), 13-58 (P13- 
P58), 59-64 (PT3-PT8), 65-96 (F1-F32). The time span of 
the diagram is 1 hour 30 minutes. 



FIG. 8: (Color online) A typical SSD for a non-focal electrode 
[here, AT5 (i=5)]. Electrodes are numbered as in Fig. [7] The 
time span of the diagram is 1 hour 30 minutes. 



represents a good compromise to obtain maximum detail 
with minimal noise. The calculations were verified with 
time windows of 1 and 20 seconds. 

The calculation is repeated for all signals as a refer- 
ence in turn, since Sij ^ Sji in general. The study of 
asymmetry is beyond the scope of this paper. 

We refer to the diagram in which Sij(t) is displayed 
as the Strength of Synchronization Diagram (SSD). We 
color code the values of Sij (t) such that black is zero and 
white is the highest. 

Two such typical diagrams [using as reference elec- 
trodes the specific electrodes F14 (i = 78) and AT5 
(i = 5)] and j from 1 to 96 are presented in Figs. [71 and [U 
The diagrams correspond to a part of the record one and 
a half hour long, containing two subsequent seizures. The 
electrodes in these figures are numbered and correspond 
to: 1-8 (AT1-AT8), 9-12 (ST1-ST4), 13-58 (P13-P58), 59- 
64 (PT3-PT8), 65-96 (F1-F32) (cf. Fig©. 

Note that a "strength of synchronization" of, say, 0.2, 
means that the signals are synchronized 20% of the time 
within the given time window. 

The SSD provides a clear scenario of synchronization 
dynamics between these two seizures. The diagrams start 
with the final stages of the first seizure (ending at about 
t ~ 20 seconds) with a high synchronization strength in 
almost all channels. Then the synchronization is lost for 
about 15 minutes. It then reappears at t ~ 1100 seconds. 
In Fig. [71 showing SSD for the reference electrode F14, 
the most strong synchronization to the group of neigh- 
boring electrodes F6-8 (i = 70-72), F15-16 (i = 79-80) 
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FIG. 9: (Color online) A typical mean pair- wise synchroniza- 
tion matrix Sij, computed from the same Sij as shown in 
Figs. [7]and[8l for the time interval (2000,2600) sec. 



and F22-23 (i = 86 — 87) is evident, while a larger group 
of electrodes is synchronized with F14 at lower strength. 
Interestingly, this electrode demonstrates similar level of 
synchronization strength with a spatially remote P grid, 
especially with electrodes P49-PG52 (i = 49 - 52). At 
the same time, the electrode AT5 (FigJSJ) becomes weakly 
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FIG. 10: (Color online) The geometric spread of the elec- 
trodes, for the patient A. The viewing angle was chosen to 
best present the resected area. Only relevant electrodes are 
labeled. The electrodes with synchronization strength exceed- 
ing the threshold value S are marked by a rectangle. The 
electrodes chosen for resection are highlighted. 



synchronized with only few neighboring electrodes from 
the same strip. The next seizure occur between t « 3000 
seconds and t w 3240 seconds. Starting around t — 2940 
with a decrease in synchronization the event proceeds to 
exhibit high synchronization between all the electrodes, 
starting around t = 3000. After t « 3240 the syn- 
chronization is again lost and regained by F14 at about 
t = 4300 seconds, as it did at t « 1100. 

It is crucial to stress that the SSD patterns presented in 
Figs. and [8] are qualitatively preserved over all seizures 
of the same patient in our data. Moreover, the groups of 
strongly synchronized electrodes are preserved through- 
out all records, although at varying strength. The desyn- 
chronization following the seizure, apparent in Figs. [71 
and El corresponds to the post-ictal state, detected at 
the time of recording. 

Similar behavior was found in the data of two other 
patients. In a patient A's data only four pairs of elec- 
trodes remain synchronized during all available records. 
Two of these pairs were placed at the temporal lobe. 
During the seizures, a particularly strong synchroniza- 
tion was found in the signals coming only from AT and 
ST electrodes. For patient C three independent groups 
of strongly synchronized electrodes were found, while the 
overall behavior was similar to that of patient B. 



III. RESULTS AND DISCUSSION 

Previous observations involving different linear and 
non-linear measures [l3|, [3, HH, [22, 23] indicate that the 
values of the corresponding measures are altered over 
tumors and regions thought to be part of the epilep- 
togenic zone. The regions of such a locally enhanced 



FIG. 11: (Color online) The electrodes scheme for the patient 
C. The viewing angle was chosen to best present the resected 
areas. Only relevant electrodes are labeled for clarity. The 
electrodes are marked as in FigfTOl 

or suppressed activity are stable and have well defined 
edges. Since we found high synchronization strength 
during epileptic seizures, we identify the foci of epilep- 
tic activity with those regions that statistically preserve 
a high degree of synchronization at all times. Thus our 
definition of the foci of epileptic seizures is the ensemble 
of those electrodes that, besides being strongly synchro- 
nized to each other during the ictal period, remain also 
synchronized above a chosen threshold during the inter- 
ictal period. To identify this ensemble, we now compute 
the time average of Sij according to 

T 

Sij = T" 1 ^2 Sij(tk) , (3) 
t k =i 

where the time interval (1,T) does not include a seizure. 
An example of the resulting mean pair- wise synchroniza- 
tion matrix Sij is shown in Fig. [9j Si j converges well for 
T > 60 seconds. Clearly, the synchronization strength 
between electrodes within a given grid is higher than be- 
tween spatially remote electrodes. For a collection E re f 
of iV r ef reference electrodes and a collection £ e \ of N e \ 
electrodes (cf. Fig. [9]) we define the mean S as: 

For the diagram shown in Fig. [9l we calculated S and 
the standard deviation <j of Sij for each of five grids and 
between the grids. The value of S within the grid, av- 
eraged over five grids is 0.056 (varying 0.047 — 0.065 
for different grids). The standard deviation a within 
the grid, averaged over five grids is a — 0.034 (rang- 
ing 0.03 - 0.042 for different grids). The value of S 
averaged over all inter-grid combination of electrodes is 
0.038 (range 0.031 - 0.045) within a = 0.007 (range 
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0.004-0.0125). At the same time, the SS of the strongly 
synchronized group ranges in the interval 0.1 — 0.2, be- 
ing more than three standard deviations above the corre- 
sponding mean value. The difference is especially strik- 
ing for the inter-grid synchronization, found between the 
group in F grid (F6-7, F14-16), and the remote group in 
P grid (P49, P50). These values are typical for all data 
of patient B, except for post-ictal periods. Similar values 
of Sij were found for patients A and C, although no such 
remote synchronization was detected. This observation 
is consistent with the relations between values of mean 
phase coherence of the electrodes pairs that belong to the 
local hypersynchrony regions and those that are outside 
of these regions, found in [14j . 

To localize the foci we need to choose a threshold value 
5th, such that only electrode pairs with Sij > 5th will be 
selected as belonging to the focus. Although for these 
pairs usually both Sij and Sji were above the threshold, 
we always used the largest of the two to compare with 
5 t h- Clearly, the number of selected electrodes is a strong 
function of 5 t h- To compare with the actual selection we 
chose 5th, such that the number of electrodes designated 
by the present procedure agrees as much as possible with 
the number chosen by the epileptologists. It turned out 
that chosen values of 5 t h are larger than the value of 5 
by at least three standard deviations. 

The results, that refer to calculations over the entire 
set of available data, are displayed in Figs. [2j [10] and 
fTTI We color those electrodes that were chosen for the 
operation, and we put a rectangle on every electrode that 
is identified by our method. 

For patient A, the SS was relatively low in the inter- 
ictal periods. On the other hand, during the seizure, the 
entire region, covered by AT and ST strips was synchro- 
nized at a level of 0.6-0.7. We included in the selection 9 
electrodes, that showed highest pairwise synchronization. 
They all are placed within the resected area, see Fig. [10J 
This patient was seizure-free after the surgery. 

We have shown in Figs 09] that patient B had two 
remote groups of electrodes that remained inter-ictally 
synchronized. As seen in Fig [2j the chosen subset of 
strongly synchronized electrodes overlap with both foci, 
resected during surgery. Out of 13 electrodes marked 
for resection, 9 were identified also by us. However, we 
have found also other electrodes, that belong to the same 
focus (in the frontal grid), as well as adjacent electrodes 



from the parietal grid. The outcome of the surgery was 
a improvement in patient's condition, who however did 
not become seizure-free. 

As for patient C (cf. FigfTT]). the operated area and 
our chosen set of electrodes overlap only partially. Out 
of 30 electrodes marked by epileptologists, only 9 were 
chosen also by us. An additional strongly synchronized 
group was found at the edge of the parietal grid, indicat- 
ing the possibility for another, not resected focus. After 
initial improvement, this patient did not benefit from the 
surgery. 

In summary, we pointed out that a simple and com- 
putationally fast method to detect phase synchrony in 
ECoG signals was able to detect and localize the mul- 
tiple foci of epilepsy. The method does not require any 
assumptions about the spectral features of the signal or 
the nature of the synchronization between the signals. 
The measurement of the synchronization strength from 
the synchrogram amounts to a straightforward count of 
the data points within a chosen interval. The synchro- 
nization scenario in the long recordings is characterized 
by repeated patterns from one seizure to the next one. 
For each patient, we found a subset of electrodes, that 
statistically preserve a high degree of synchronization at 
all times. These subsets overlap fully or partially with 
the areas chosen for resection. In patient B we correctly 
identified the existence of two distant foci. In patient C 
we also identified the presence of additional focus, not 
operated on. The overlap of the location of the elec- 
trodes identified by our method and the resected areas 
was consistent with the degree of surgery success. 

Although it is not reasonable to draw any final con- 
clusions on the basis of the study of only three patients, 
these observations suggest that this method may become 
a useful tool to help improve the outcome of the surgical 
treatment of epilepsy. 
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